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Abstract 

We present the first complete calculation of the one-loop electroweak effect in the process of semi- 
. 1^ inclusive bottom-Higgs production at LHC in the MSSM. The size of the electroweak contribution 

X 

;_i depends on the choice of the final produced neutral Higgs boson, and can be relevant, in some 

range of the input parameters. A comparison of the one-loop results obtained in two different 
renormalization schemes is also performed, showing a very good NLO scheme independence. We 
further comment on two possible, simpler, approximations of the full NLO result, and on their 
reliabilty. 

PACS numbers: 



I. INTRODUCTION 

It is a well known fact that tan P enhanced Yukawa coupling in the Minimal Supersym- 
metryc Standard Model (MSSM) could favour the Higgs production in association with 
bottom quarks, contrarily to the Standard Model (SM) case, where the Higgs production is 
dominated by top-Higgs coupling. 

Due to its relevance as a possible channel for the Higgs discovery, in the last few years 
the associated bottom-Higgs production has been extensively studied in the literature. 
Depending on the choice of the flavour- scheme in the partonic description of the initial 
state and on the identified final state, one can consider a number of different partonic 
sub-processes for Ti^ + bjets production: while the choice of the 4 versus 5 flavour scheme is 
mainly theoretically motivated, resulting in a reordering of the perturbative expansion [1] , 
the requirement of a minimum number of tagged b in the final state is physically relevant in 
the signal extraction. Assuming the 5-fiavour scheme (which ensures a better convergence 
of the perturbative series resumming large logarithms in the bottom PDF), one can consider 
three different types of production processes, depending on the required final states: the 
exclusive one where both bottom jets are tagged {hlfHP final state), the semi-inclusive 
one where only one bottom quark is tagged {hl-L^), and the inclusive one where no 
bottom quark jets are tagged. While the inclusive process has a larger cross section [2], the 
semi-inclusive with a high p^ ^ bottom in the final state is experimentally more appealing [3] . 

The relative weights of the partonic processes {bb — > "H^, hg — )■ 6?^°, gg — > WH^) are 
analyzed in [2], where also the Ug corrections (NLO) to the leading sub-process hb — t- 'HP 
are computed. The NNLO order in QCD (a^) for the same sub-process is calculated 
in |3], while the electroweak (SM and MSSM) and SUSY-QCD NLO corrections have been 
computed in [S], showing that the size of electroweak corrections can be comparable, for 
large tan/3, with that of the strong ones. 

The associated semi-inclusive production process {IfHP final state) is analyzed at the 
NLO in QCD in [3J and |iij, while the effect of the SUSY QCD is given in [7j. Very recently, 
Dawson and Jaiswal have also computed, for the Standard Model process hg — )■ bhsM, the 



one-loop weak corrections [8]. 

Finally, the exclusive process, where two bottom jets are tagged in the final state, is 
considered at the NLO in QCD in [1], [9], ^U\ and [H]. The leading Yukawa corrections for 
this partonic process are considered in [1^ and SUSY QCD effects have also been computed 

in na. 



Our paper is strongly motivated by the possible relevance of the associated bottom-Higgs 
production in the experimental search of the Higgs at the LHC; moreover, as stressed in [3], 
the susy one-loop ew effects (for the inclusive process) can be sizable and they can be safely 
accounted by an improved born approximation. Therefore the spirit of our computation is 
twofold: on the one hand we provide for the first time the complete NLO EW corrections 
for the semi-inclusive process, including also the overall QED effect, that was not computed 
by [H] , and on the other hand we can perform a further and independent test on the validity 
and limits of the improved born approximation in different scenarios. Our calculations have 
been performed in two different (DR and DCPR) renormalization schemes: as expected the 
final one-loop results are, within at most a relative few percent difference, the same in the 
two frames; however, the DR scheme appears to be the one where the perturbative effect is 
numerically mostly more under control. Therefore we shall discuss our results in this frame, 
showing in various figures the dependence of the different observables on the choice of the 
input parameters. We have finally compared the results obtained with the full electroweak 
computation with those obtained within a commonly used approximation scheme. This will 
be done in the final part of our paper, which is organized as follows: Section [IT] contains 
a general concentrated discussion of the actual derivation of the theoretical formulae (a 
part of which has been shifted in a technical Appendix |B|) to be used for the calculation of 



the various observables. Section III and IV contains our numerical results, that are briefiy 
discussed in the final Section |Vl 



II. KINEMATICS AND AMPLITUDE OF THE PROCESS bg -^ bVP 

A. Kinematics 

At lowest order there is only one partonicP channel leading to bottom-Higgs production 

KPb) 9{Pg) ^ KK) ^°(Pwo) (1) 

where H^ is one of the three MSSM neutral Higgs bosons (/i°, if°, A^). In the partonic 
center of mass frame the momenta of the particles read 

Pfc = (^fe;0,0,p), pg = (p;0,0, -p), 
p'fe = (E^; p sin e, 0, p cos 9), p^o = (E^o; -p sin 9, 0, -p cos 6'). (2) 

The Mandelstam variables are defined as 

s = {pg+ Pbf, t={pb- p'J^ u = {pg- p'f,y. (3) 

For later convenience we define two momenta q and q' as follows 

<l = Pb+P9, q' = p'b-Pg 

B. Born and one-loop amplitudes 

We denote the 0{a'^a'^) contribution to the amplitude (differential cross section) of the 
process X as Ai'^ {da'^). The Born terms result from the s- and u-channel h quark exchange 
of Figure [TJ The color stripped tree-level amplitude reads as follows 

-Mjg^ili'o = -f^H)4(^6)[c'^(&&^°)^L + c^(&&^°)Pij](^ + m,)/j;/x)«,(A,), 

^"j 4(A'J/j:/i)(/+ mb)[c\hh'H')PL + c\hhH')Pn]ub{\b), (4) 



u — nif^ 



where A;,, (A'f,) is the helicity of the initial (final) bottom quark while /i is the polar- 
ization of the gluon. Ub{\b) WbiK)] i^ the spinor of the initial [final] bottom quark, 



[1] One should also consider the photon induced process bj — > bH'~': the contribution to the total cross section 
arising from this sub-proceess is doubly suppressed, due to the smaller 7 parton distribution function and 
smaller coupling {a instead of as). Since the main goal of this paper is the calculation of the NLO 
electroweak effects for bH'^ production, and the bj —>■ bT-L^ can be safetly computed at the LO, we do not 
take into account the photon induced production in the following. 



ec,(/i) = (0; /i/v^, — i/v^, 0) is the gluon polarization vector and Pr,l = (1 ± 7^)/2 are 
the chirahty projectors. The relevant couplings c'^{bb'H^) [rj = L,R) are defined as 

errih \ sin a 



^ ' \2swMwJ cos/3' ^ ^ V 



\2swMw J cos/3 



c^(66A°) = -t (-^^^] tan/3, c^(6M°) = c^*(6M°) . (5) 

\2swM\Y J 

We factorize out of the gluon couplings the colour matrix element A°/2. The sum over colors 

leads to a factor 

\^ fA'^ AM , 



.2 2 , 

a=l ^ ^ 

that multiphes the squared amplitude. 

The generic helicity amplitude can be decomposed on a set of eight forms factors J''^ {rj = 

L,R) as follows 

\k=lri=L,R ) 

where 

The only non-zero scalar functions at the tree level are -^fj^^^-^o and A^^\.5->^o- They read as 
follows 

og^DH S — mf^ U — ml og^OH U — TTl^ 

The one-loop electroweak virtual contributions arise from self energy, vertex and box 
diagrams. Counterterms for the various bottom quarks lines, for the "H" line, and for the 
bbTi^ coupling constants have to be considered as well. The corresponding diagrams can be 
read off from Fig |2j |3| 

All these contributions have been computed using the usual decomposition in terms of 
Passarino-Veltman functions and the complete amplitude has been implemented in a C-|--f 
numerical code. 



C. Renormalization 

In order to cancel the ultraviolet (UV) divergences the Higgs sector and the bottom 
sector have to be renormalized at 0(a). The expressions of the counterterms entering our 



calculation are collected in Appendix [Bj 

Higgs sector 

As anticipated we performed the calculation using two different renornialization schemes: 
the DR scheme [l4j is defined by the following renornialization conditions 



<5Z- 



^^h: 



Re 



9E^o(A;2 



fc2=M2.,a=0 



9P "'=''"hO 






STho = -Tho 

SMlo = Re^AoiMlo) - MloJ:'^o{Mlo) 
6 tan /3^ = - {5Z^^ - 5Z^) tan /3. (10) 

SZ^ define the wave function renornialization costant of the Higgs field Hi, the third and 
fourth line fix the tadpole renornialization and the last one the tan/3 renornialization con- 
stant. [^]div means keeping the UV divergent part of A, discarding the finite contribution. 
In the DCPR scheme fi5[ [16] the independent parameters are the same, and the renornial- 
ization conditions of the Higgs wavefunctions change as follows 



(9P i«=^%o tan (3 Mz 

.^DCPH _ p /S^o(fe') | , tan/3 2^ 

dZjj^ - -Re — —^ — \k^=Mlo + -JJ^^^^AoziM^o) 

6ThO = -T/jO 

6Tho = -Tho 

SMlo = Re^AoiMlo) - M%X'ao{MIo) 

S tan /3°^™ = ^ (5Z°f « - 5Z^f '') tan /3 (11) 

We choose to impose on-shell (OS) condition for the mass of CP-odd A^ Higgs in both 
schemes. 

The renornialization constants of the Higgs bosons wavefunctions and of the d'{hhl-i^) 
couplings can be written in terms of the of the renornialization constants defined above. 



Their explicit expression is given in Appendix \K\ 



Bottom sector 

The mass of the bottom and its wavefunction renormahzation function is fixed in the on-shell 
scheme: 

(5m°^ = ^mfe [ReEfeJm^) + ReE^^K) + 2ReSb^K)] , (12) 



6Zl; = -ReSb.K) - ^'^Re [S.Jfc^) + ^^^(A;^) + 2Sb,(P 
6Z^ = -ReSbJm2)-m2ARe[s,,(A;2) + S,,(A;2) + 2S,,(A;2 



Ifc2=m? ! 



lA:2=mj 5 



where the bottom self energies are defined according to following Lorentz decomposition: 

Mp) = ^PL^bAP^) + l^PR^bAp") + rn,T.,Ap') ■ (13) 

The bottom masses in the Yukawa couplings are treated completely at the electroweak level, 
with OS or DR renormahzation conditions respectively in the two schemes. Resummation 
of large logarithms from the running of the bottom mass suggests to trade bottom mass 
appearing in the couplings with an effective bottom mass, |1TJ . The resummation of the 
(astan/3)" contributions can be achieved modifying the tree level relation between the bot- 
tom Yukawa coupling and the bottom mass: the bottom mass of the couplings, which is 
related to the bottom Yukawa coupling, is replaced by an effective mass, (e. g. in the DR 
scheme) 



mT -> rnT = ^f^ (14) 



where A;, is given by 



2 a 

A, = -^ Mgf, tan 13 I{M-,^,M-,^,Mg) (15) 

on 

Moreover, the bbHi coupling is dynamically generated at 0{as) and can be enhanced if 
tan/3 is large. This effect can be included modifying the c^{hh'HP) couplings. The actual 



effect of this modification and of the bottom mass resummation, Eq: (14), is to substitute 



the c^{hh'hP) couplings in Eq. ^ as follows 



n, . C\hhh^) ^^ _^ [^ A, 



d^ihhK') -> ^ '- X m°" 1 - 

mfc \ tan p tan a 

v^i.r.uO^ c"{bbH^) _^ / , Afctana 
d'ibbH') -)■ — ^ X m " 1 + 



r/if, V tan /3 



c''(6M'^) -^ 



c 



(^x^-fl-^V (16) 



We have checked the cancellation of the UV divergences among counterterms, self-energies 
and triangles. This cancellation occurs separately inside 8 sectors, i.e. s-channel "initial" 
triangles with chirality L or R, s-channel "final" L or R, w-channel up triangles (L or R) 
and M-channel down triangles (L or R). The Box diagrams are UV finite. 

D. QED radiation 

The infrared (IR) singularities affecting the virtual contributions are cancelled including 
the bremsstrahlung of real photons at Olaso"^), 

KPb) g{Pg) -> h{p',) n^pno) lip,) , (17) 

arising from the diagrams in Figure |4j This contribution has been computed using 
FeynArts ^18j and FormCalc ^19j. The integral over the photon pase space is IR divergent 
in the soft-photon region, i.e. for p° — )■ 0. The IR divergences are regularized within 
mass regularization, giving a small mass m, to the photon. The phase space integration 
has been performed using the phase space slicing method. This method introduces a 
fictitious separator AE and restricts the numerical phase space integration in the region 
characterized by p^ > AE. The integral over the region p^ < AE is performed analytically 
in the eikonal approximation 



Large collinear logarithms containing the bottom mass can be re-absorbed into the redef- 
inition of the parton distribution function (PDF) of the bottom fb{x,fi). In the MS (DIS) 
factorization scheme this is achieved performing the substitution [21] 



fb{x,fi) -^ fb{x,fi) <{ 1 el 



1 - In 5s - In 5, + ( In (5s + - ) In ( — ^ ) - ^Apc/^i 



i)'°(5)-r- 



a 2 [^ " dz „ /x 
2^""' L ~z 



fb(-,^^ 

\z 



l + z\ fl' 1 \ 1 + z' , 



8 



,(1^ 



and setting Apc = (Apc = !)• /i is the factorization scale, Ss = 2AE/y/s, while Cb is the 
bottom charge, ni and k,2 are defined as follows, 

2 



1 + ^2 /i-z\ 3 1 



K2 



I- Z 



In 



21 -;z 



+ 2z + 3. 



(19) 



We tested numerically the cancellation of IR divergences, the independence of our results of 
m^ (in the sum of the soft and virtual part) and of the separator AE (see Figures [sl [gI [t]). 



E. Total cross sections 

Including the finite wave function renormalization for the Higgs field we obtain the fol- 
lowing expressions for the tree-level differential partonic cross section of the processes we 
are considering, 

P'dcOsO ^ , A/2,1/2 2 , 



da 



1.1 



bg^bno 768 7rs/3 
where (3 = 2p/y/s, (3' = 2p' /y/s, and s is the 



Z-HO 



A.l/2,1/2 



Mandelstam variable defined in Eq. pi); the 
NLO-EW contribution to the differential cross section reads as follows 



,.1,2 _ f3'dcos9 



i-z. 



cos a 



hOHO 



sma 



_A^l/2,l/2 
■'^^bg^bhO 



+ 2ReAi;/!;'£.(A<;«^.)-}-d*;;_ 



•bftO' 



da 



1,2 



f3' d cos 9 
^3^^"' ^ 76871SP 



'HO- 



1-Z 



sma 



HOhO ■ 



COS a 



_A^l/2,l/2 









(3' d cos 
7687rs/3 



(21) 



we 



where the Z factors Z/^o, Z^o, Z^o, Zf^of^o, and Zfjo^o in the two renormalization schemes 
are considering can be found in [H] and in [TJ]. The partonic differential cross section for 
the real photon radiation process reads as follows. 



,.1,2 i 

"%^6W07 "4-24 



d<t>{p'b,PHO,Pj) 
2/3 s 



^■H" 



■^^^g^bU^-y 



(22) 



where, according to the notation introduced in [22], (i0(p'^,p-^o,p^) is the three-particles 
phase space measure. The hadronic differential cross section at 0{asa) and (^(a^Q;^) reads 



da^p^^^y^a{S) 



7 1,2 



(S) 



^1 .1 

dxi / dx2 
'0 Jo 

.1 .1 

(ixi / (ix2 
'0 Jo 



X 



.1,2 



fb{xi,fi)fg{x2,fi) + (Xi ^ X2) 



fb{xi,fi)fg{x2,fi) + (Xi 4-> X2) 



1,2 



c?^b„^fea.o(a;ia;25) + da'.^oJxiX2S) 



' bg^bW 



bW-y^ 



(23) 



respectively. yS is the hadronic center-of-mass energy, while fi{xi, jj) is the parton distri- 
bution function of the parton i inside the proton with a momentum fraction Xj at the scale 
/i. For later convenience we define the invariant mass distribution as 
-1,1 ^1 ri 



da 



pp^bn° 



^1,2 
"'"pp^bu^ 

d^/% 



dxi / dx2 
Jo 



fb{Xi,fi)fg{x2,fi) + (Xi <->X2) 



X da 



1,1 



o(xiX2S')d^ ( a/xiX25' 



X 



bg^bU 

1 .1 

(ixi / (ix2 
io 

1,2 



/fe(Xi,/i)/g(x2,/i) + (Xi ^X2) 



da 



bg^bn 



o{XiX2S) + C?6-^^^^_^^^o^(XlX25') 6 \^^/xiX2S - y/sj 



(24) 



III. NUMERICAL RESULTS 

The independent input parameters in the MSSM Higgs sector are the A" mass and tan/3: 
since we impose the same renormalization condition for M^o only tan /3 should be converted 
in the change of scheme, using the one-loop relation: 



tan/3^ 



tan /3°''- + 5 tan /3°" - 5 tan /3^ 



(25) 



while the OS and DR bottom masses m^^ and m™{fi) are computed starting from 
m^^(mf,) = 4.2 GeV and following the procedure described in Section 3.2.2 of [T7]. 

For the numerical evaluations we used the supersymmetric scenario SPPi and a class 
of points of the parameter space SPP2, with variable tan/3 = 10,20,30,40. The input 
parameters characterizing these scenarios are summarized in Table [Tj The sparticle masses 
and mixing angles have been obtained with the code FeynHiggs [23]. The one-loop Higgs 
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Scenario 


tan/3 


M^o 


Mq,i 


Mg,2 


^9,3 


Ml 


M2 


Mg 


SPPi 
SPP2 


15 

variable 


350 
250 


350 
500 


350 
500 


250 
400 


90 
90 


150 
200 


800 
800 



TABLE I: Inputs parameters for the SUSY scenarios considered in our numerical discussion. Mqj 
is the common value of the breaking parameters in the sector of the squarks belonging to the j"' 
generation. The dimensionful parameters are given in GeV. 



tan^ 


DR ,NLO 


_DR ,LO 


DCPR,NLO 


DCPR.LO 


^DR 


J^DCPK 


NLO ratio 


10 


1.367 


1.281 


1.371 


1.253 


1.067 


1.093 


0.997 


20 


5.040 


4.784 


5.060 


4.278 


1.053 


1.182 


0.995 


30 


10.601 


10.295 


10.785 


8.505 


1.029 


1.268 


0.98 


40 


17.118 


17.125 


17.615 


13.038 


0.999 


1.350 


0.97 



TABLE II: A^ production, SPP2 spectra: total cross sections [pb], i^- factors and NLO DR /DCPR 
ratio 



tan/3 


_DR ,NLO 


^DR ,LO 


_DCPR,NLO 


_DCPR,LO 


^DR 


^DCPR 


NLO ratio 


10 


1.338 


1.260 


1.340 


1.234 


1.061 


1.086 


0.998 


20 


5.133 


4.857 


5.099 


4.334 


1.056 


1.176 


1.006 


30 


10.975 


10.461 


10.715 


8.488 


1.049 


1.262 


1.024 


40 


18.613 


17.918 


17.811 


13.248 


1.038 


1.344 


1.045 



TABLE III: H^ production SPP2 spectra: total cross sections [pb], K-factors and NLO DR /DCPR 
ratio 



tan/3 


DR ,NLO 


DR ,LO 


DCPR.NLO 


DCPR.LO 


^DR 


TV-DCPR 


NLO ratio 


10 


0.282 


0.248 


0.282 


0.243 


1.135 


1.156 


1.002 


20 


0.255 


0.254 


0.254 


0.230 


1.005 


1.107 


1.003 


30 


0.228 


0.258 


0.230 


0.217 


0.882 


1.059 


0.988 


40 


0.204 


0.267 


0.213 


0.211 


0.764 


1.012 


0.955 



TABLE IV: h^ production SPP2 spectra: total cross sections [pb], i^-factors and NLO DR /DCPR 
ratio 
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masses are numerically computed by finding the zero of inverse one-loop propagator matrix 
determinant 

[e - Ml, + S^o ( A;2)] [k^ - Ml, + S^o ( A;2)] - S^o^o ( A;^) = 0. (26) 

Since we require semi-inclusive production (i.e. the bottom quark must be tagged) we 
impose the following kinematical cuts on the bottom in the final state, limiting the trans- 
ferred momentum pb^^ > 20 GeV (due to resolution limitations of the hadronic calorimeter) 
and the rapidity \yh\ < 2 (in order to be able to perform inner tracking). The process we 
are considering is leading order in QCD. Therefore, analogously to [21H26] . we use a LO 
QCD PDF set, namely the LO CTEQ6L [27J. Our choice is justified since the QED effects 
in the DGLAP evolution equations are known to be small [2H]. The factorization of the 
bottom PDF has been performed in the DIS scheme, with factorization scale /i = M-^o + m^^. 

In Figures [SJ [9, 10 we show the total cross section for A^, if° and h^ production in the 
class of supersymmetric scenarios SPP2, as functions of tan/3. We present both the results 
in the DR and in the DCPR schemes. The numerical values and the K-iactois in the two 
schemes (defined as usual as the ratios a'^^^ /a^^\ note that the LO is computed with the 
resummed/modified SUSY QCD coupling, so our i^-factors account of the pure electroweak 
NLO effect), as well as the ratios of the NLO cross sections in the two scheme are reported 
in Table [TTl [TTll |IV 



As one sees, the values of the total cross sections do coincide in the overall range, apart 
from small differences of the few percent size for very large tan /3 values. This confirms our 
expectation that at the NLO level the two schemes should be equivalent, and also provides 
an important check of the reliability of our calculations. 

Having verified the realistic one-loop equivalence of the two schemes, we have decided 
to perform our analysis in the DR scheme. The main theoretical reasons of our choice 
have been fully illustrated in [22] • In particular this scheme is known to be generally more 
stable numerically: our results confirm mainly this expectation but it is worth to note that 
for h^ production both schemes can produce (in different tan/3 regions) relatively large 
effects; nevertheless the good agreement between the two schemes leads to suppose that the 
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H" 


_DR ,NLO 


^DR ,LO 


^DR 


^0 


0.768 


0.724 


1.060 


H^ 


0.769 


0.727 


1.056 


h^ 


0.213 


0.222 


0.961 



TABLE V: SPPi spectrum: total cross sections [pb] for the three Higgs and DR if-factors 
perturbative expansion is well behaved, and NNLO effects are well under control. 



Figure 14 shows the i^-factors for the three Higgs bosons in DR as function of tan /3 while 



Figures 15, 16, 17 show, for the scenario SPP2 tan/3 = 30, the invariant mass distribution 



and the relative NLO effect. In the next Figures 18, 19, 20 we again plot the differential 



distributions for the SPPi scenario; the total cross sections for this scenario are reported in 
Table IVl 

From inspection of the figures, one can draw the following main conclusions: 

1. The i^-factors for H^,A^ are systematically small for large tan/3, and would reach a 
larger size (roughly, 8 %) for small tan/3 values around 10. 

2. The K-factoT for h^ varies drastically with tan/3, changing from positive values of 
about 15 % for tan/3 around 10 to negative values of about 25 % for tan/3 around 40. 
These extreme negative and positive values are of a size that cannot be ignored in a 
dedicated experimental analysis. 

These features follow from the THDM structure and the h^ — H^ mixing where a is close 
to /3 — tt/2 leading to a tan/? enhancement in the h^ case but to a l/tan/3 suppression in 
the if °, y4° cases. 

This, we believe, is the main message of our calculation: while for sure the QCD NLO 
are the dominant corrections (of order 20 — 40% depending on the Higgs mass, see for 
example |3j), as it was to be expected from the analysis of Dittmaier et al. [5], the one-loop 
electroweak contribution in the semi-inclusive bottom-Higgs production processes must not 
be a priori considered as negligible. 
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IV. NUMERICAL APPROXIMATIONS 

Having performed the calculation of complete one-loop effect on the process, we shall 
consider the possibility of simpler, effective approximations to the full and long calculation, 
that may be used to obtain a quicker and qualitative description of the results. 

With this purpose we have first considered the "improved Born Approximation" (IBA) 
following the prescriptions given in [5j: the IBA is obtained is this case by including in the 



definition of A^ (see eq. 14) the electroweak contributions and then replacing the mixing 
angle a with the effective value «£//, obtained by the diagonalization of the one loop mass 
matrix 

The effect of the latter redefinition of a is negligible for H^ and A^, but significant for h^. 



As one can see from the plots, (Figures [lT|l2|13[ ) this version of IBA is sufficiently close 
to the complete calculation only for relatively small tan/3 values, roughly tan/3 < 20. In 
this range, the approximation gives larger (compared to the complete calculation) rates for 
H^, A^ and smaller rates for /i°. The differences remain below the ten percent size, which 
would be tolerable at least in a first phase of LHC measurements. Increasing the tan/3 
value, the IBA description becomes worse. For tan/3 = 40, it differs in all the three cases 
by, roughly, a relative 25 percent, which seems a rather poor prediction for the measurable 
total rates. 

For what concerns the tan/3 dependence of the plots, one can conclude that it provides 
those features that would be expected at the chosen value of M^o, which is sufficiently 
larger than Mz to approach the correct decoupling limits. In this large M^^o regime, that 
is discussed widely in the literature (see e.g. [SO]), the iJ° and A^ couplings become almost 
exactly proportional to tan/3, while the h^ coupling becomes very weakly tan/3 dependent. 
These features are well reproduced by the plots, that show a roughly quadratic tan/3 
dependence of the if°, A^ rates and a much weaker tan/3 dependence for /i°. But for large 
tan /3 values, there seems to be an extra tan /3 dependence of the complete calculation that 
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tan/3 


^DR ,NLO 


RVAhbn 


a/RVAbbn 


IBA 


a/lBA 


10 


1.338 


1.32623 


1.00888 


1.34087 


0.997861 


20 


5.133 


5.08324 


1.00979 


5.48397 


0.936 


30 


10.975 


10.8433 


1.01215 


12.6044 


0.87073 


40 


18.613 


18.3461 


1.01455 


22.6229 


0.822749 



TABLE VI: H^ production: comparison between the complete NLO prediction and the two ap- 
proximations: total cross sections and ratios u^'^ '^^°/cJapp.. 

is not contained in the IBA description. 



Having this apparent discrepancy in our mind, as a second attempt, we have tried to 
use what we would call a "Reduced Vertex Approximation" (RVA): we approximate the 
complete NLO keeping only the (all) one loop corrections to the "final" Yukawa bb'H'^ 
vertex and the subset of counterterms needed to get a UV- finite result; the photon mass is 
regulated (arbitrarily) as M^ = Mz (and thus we do not include soft and hard radiation). 
We kept the one loop Higgs masses in the kinematics as well as the Z-factors in the 
definition of the cross section; all the other diagrams (Boxes, Initial and Up Triangles, Self 
Energies) are neglected. As a check we computed the cross section in this approximation 
in both schemes (the subset of diagrams, with the right choice of counterterms, should 
be scheme independent). As one can see from the updated figures our RVA turns out to 
provide very efficient description of the total NLO cross sections; the difference between 
the NLO and the RVA is of order of 1%, 3.4% in the worst case. This is numerically 
summarized in Tables |VT]|VlT]|Vni] and Figures [TTp2p3 



From the inspection of those Tables and Figures we would conclude that the extra vertices 
that the RVA contains seem to provide the extra tan (5 dependence not predicted by the IBA 
in a reasonably satisfactory way, i.e. at the level of few percent in the full tan /3 range. This 
RVA cannot be transformed into simple analytical expressions. It tells us that the relative 
effect of a large set of Feynman diagrams, those that were not included in the approximation, 
is small, at the level of a few percent, which might be considered negligible for the first phase 
of LHC measurements. 
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tan/3 


_DR ,NLO 


RVA,,^ 


a/RVAbb^ 


IBA 


cj/IBA 


10 


0.282 


0.277157 


1.01747 


0.268161 


1.05161 


20 


0.255 


0.250495 


1.01799 


0.238459 


1.06936 


30 


0.228 


0.221673 


1.02854 


0.211275 


1.07916 


40 


0.204 


0.197159 


1.0347 


0.164874 


1.23731 



TABLE VII: h^ production: comparison between the complete NLO prediction and the two ap- 
proximations: total cross sections and ratios 0"°^ ■'^^^/cjapp.. 



tan/3 


_DR ,NLO 


RVA56^ 


a/RYAbbn 


IBA 


ct/IBA 


10 


1.367 


1.35328 


1.01014 


1.36737 


0.999729 


20 


5.04 


4.98026 


1.01199 


5.4543 


0.924042 


30 


10.601 


10.4581 


1.01366 


12.2948 


0.862232 


40 


17.118 


16.8292 


1.01716 


21.7326 


0.787663 



TABLE VIII: A^ production: comparison between the complete NLO prediction and the two 
approximations: total cross sections and ratios a°^ '^^° /aApp.- 

V. CONCLUSIONS 

We have performed in this paper a complete MSSM calculation of the electroweak 
NLO effect in the processes of semi-inclusive bottom-Higgs production. Our analysis has 
been performed for two choices of the M^o input parameter and for variable values of the 
tan/3 parameter defined in the DR renormalization scheme. Although a more extended 
analysis of the parameter space would be interesting, we have found certain results that 
appear to us to be general and worth publishing. The first conclusion is that two different 
renormalization schemes appear to be practically identical at SUSY NLO as one would a 
priori expect. Working in the DR scheme, that seemed to us to be somehow preferable, we 
have found that the pure electroweak one-loop effect in the three considered production 
processes is of a size that might be relevant and therefore that this contribution cannot be 
ignored for a proper experimental analysis of the reactions. 



There could exist simpler calculations involving a smaller (but still large) number of 
diagrams, that would provide a valid numerical result. We have seen that one possible 
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Improved Born Approximation does not reproduce the correct result in a satisfactory way. 
We have also seen that another "Reduced Vertex Approximation" (which considers only the 
1-loop correction to the Yukawa bb'H'^ vertex) appears to better approximate the full NLO. 

However, if a theoretical prediction of the total cross section is requested at the percent 
level, which might be the hopefully desirable final LHC goal, our conclusion is that the 
complete one-loop calculation of the electroweak part that we have performed in this paper 
should be considered, together with the available, large, QCD corrections, as the correct 
proposal to offer to the experimental community. 

There remains a couple of relevant points to be still investigated. The first is that of 
combining this analysis with an analogous one to be performed for the process of associated 
top-charged Higgs production, for which our group has already provided a complete one- 
loop electroweak analysis [26J. The second one is that of trying to relate the DR tan/3 
parameter, which is not a measurable quantity, to a measurable tan/3 (which could be 
defined for instance by A^ -^ t^t^ decay as suggested in [29j). This would allow to draw 
plots where also the horizontal axis represents a measurable quantity. These points are, in 
our opinion, quite relevant but beyond the purposes of our analysis; work is in progress on 
these issues. 
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Appendix A: Renormalization constants in the Higgs sector 



The renormalization constants of the wavefunction of the Higgs bosons A°, /i°, H^ and of 
the Goldstone boson G^ are given by 

SZhoho = cos^ adZni + sin^ a5ZH2, SZhoho = sin a cos a{5ZH2 — SZhi), 
SZffihO = sin^ adZffi + cos^ a6ZH2, ^^a^a'^ = sin^ (3SZhi + cos^ (36Zh2, 
SZco^o = cos (3 sin (3(6 Zh2 ~ SZhi), SZ^oao = SZho^o = SZ^ogo = SZhoqo = 0. 

(Al) 

The renormahzation constants for the c^{hhh^) and for the d^{hhH^) couphngs is obtained 
differentiating the tree-level expressions in Eq. ([5|), 

5 cos (3 ' 



Sc'^ibbH'' 

Scos(3, 6M^, and 6g, reads as follows 
S cos (3 = — sin^ (3 



5g_ Srub 

g nib 

5g_ 5mh 

9 rnf, 



2M^ 



cos/3 

6 cos (3 

cos (3 



ci{bbh^), 
c'^ibbH''). 



(A2) 



Si 
9 



2 ^Stan(3 



tan/3 
ReS^(M2,), 

2 



swCwMl 



Cw 



-swM^^^^ 



S^z(O) + 



-w 



6Mj 






SL(0) 



(A3) 



with 6M^ = ReS^(M|). The c^lbbA^) couplings depends only on the angle /3. When 
computing the the renormalization constant 6c^{bbA^), one has to distinguish between the 
/3-dependent factors originated by the Hi, H2 mixing and the /3- dependent factors from the 
Hi, H2 couplings. Only the latter have to be renormalized. In particular the factor sin/3 
[l/cos/3] entering the c^lbbA^) coupling is originated from the Hi, H2 mixing [couplings], 
and thus 6c^{bbA^) reads 

9 ^b 



Sc'^ibbA^ 



6g 6mb SM^^ 6 cos (3 



2M^ 



cos/3 



c'^ibbA^ 



(A4) 



Appendix B: Contributions of the counterterms 



In this appendix we list explicitely the contributions of the counterterms writen in terms 
of the renormalization constants introduced in Section II C and in Appendix |X] The vertices 
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counterterms can be written as follows 

4 



,fc=l ri=L,R 



where J^"^ are defined in Eq. ( 8 ) while the non-zero V^ \bn° ^^^^^ 

b \ \ / ^0 

Uz'^ + ]-5Z'\ c'^ibbn') + 5c\bbH'') + \Y. ^Z*^o-Hoc\bb'H''] 



n'A.^0 = ^^\{^i^5z'i + \5z'^c\bbn^^ 

Vbl^bno = -^ (^^? - ^^3 -^(^&^°) + -^ (^4 - ^^;) c-m% (B2) 

where (r/, r/) G {(L,i?); (i?,L)} and H°,H° = h'^,H°,A°,G'^. The bottom self energy coun- 
terterm reads as follows 

The non-zero S^^^^^^o are 

+ 9s-, ^^ { u5Z^ -m^(6Zf^-2 
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Tree level diagrams 




FIG. 1: Tree level diagrams for the partonic bg — )■ bTi^ processes. 
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b quark self energy 




hO self energy 
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FIG. 2: Bottom quark self energies, higgs self energies (only the diagonal case) and internal self 



energies. 
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Triangle diagrams (part I) 



Triangle diagrams (part II) 
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FIG. 3: Triangle and box diagrams. 
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y emission 





FIG. 4: Real photon 7 emission. 
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FIG. 5: AP production: dependence of the 0{a) soft+virtual, hard, and total sum corrections on 
the separator A^ 
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FIG. 6: H^ production: dependence of the 0{a) soft+virtual, hard, and total sum corrections on 
the separator A^ 
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FIG. 7: /i*' production: dependence of the 0(a) soft+virtual, hard, and total sum corrections on 
the separator A^ 
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A production, total cross sections DR - DCPR scliemes 
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FIG. 8: Total LO and NLO cross sections in the DR and DCPR schemes, A^ production; M^o = 250 
GeV, Vh,T > 20 GeV, |yfe| < 2. 
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H production, total cross sections DR - DCPR scliemes 
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FIG. 9: Total LO and NLO cross sections in the DR and DCPR schemes, H^ production; Mj^v, 
250 GeV, P6T > 20 GeV, \yi\ < 2. 
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h production, total cross sections DR - DCPR scliemes 
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FIG. 10: Total LO and NLO cross sections in the DR and DCPR schemes, h^ production; My^o 
250 GeV, P6T > 20 GeV, \yb\ < 2. 



32 



A production, total cross sections, NLO DR and approximations 
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FIG. 11: Comparison of the total NLO cross sections: NLO DR , RVA (DR and DCPR) and 
Improved Born Approximantion (IBA), A^ production; M^o = 250 GeV, pb^T > 20 GeV, \yi,\ < 2. 
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H production, total cross sections, NLO DR and approximations 
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FIG. 12: Comparison of the total NLO cross sections: NLO DR , RVA (DR and DCPR) and 
Improved Born Approximantion (IBA), H^ production; My^o = 250 GeV, p^.T > 20 GeV, |yt| < 2. 
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h production, total cross sections, NLO DR and approximations 
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FIG. 13: Comparison of the total NLO cross sections: NLO DR , RVA (DR and DCPR) and 
Improved Born Approximantion (IBA), h^ production; My^o = 250 GeV, pb^x > 20 GeV, \yh\ < 2. 
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FIG. 14: i^-factors for A^ , H° and /i° production, DR scheme. M^o = 250 GeV, Pb,T > 20 GeV, 
IVbl < 2. 
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A production - differential distribution - DR sclieme 
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FIG. 15: Invariant mass distribution, A^ production, DR scheme. M^o = 250 GeV, tan/3 = 30 
Pb,T > 20 GeV, \yb\ < 2. 
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H production - differential distribution - DR sclieme 
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FIG. 16: Invariant mass distribution, H^ production, DR scheme. M^o = 250 GeV, tan/3 = 30 
Pb,T > 20 GeV, \yb\ < 2. 
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h production - differential distribution - DR sclieme 




0,4 0,5 

M [TeV] 

inv >- ■■ 



FIG. 17: Invariant mass distribution, /i" production, DR scheme. M^o = 250 GeV, tan/3 = 30 
Pb,T > 20 GeV, \yb\ < 2. 
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A production, differential distribution - DR sclieme 
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FIG. 18: Invariant mass distribution, A^ production, DR scheme. M^o = 350 GeV, tan/? = 15 
Pb,T > 20 GeV, \yb\ < 2. 
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H production, differential distribution - DR sclieme 
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FIG. 19: Invariant mass distribution, H^ production, DR scheme. M^o = 350 GeV, tan/3 = 15 
Pb,T > 20 GeV, \yb\ < 2. 
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h production, differential distribution - DR sclieme 
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FIG. 20: Invariant mass distribution, /i" production, DR scheme. M^o = 350 GeV, tan/3 = 15 
Pb,T > 20 GeV, \yb\ < 2. 
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